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A  thermodynamic  model  is  developed  to  determine  the  fuels  that  would  yield  an  identical  maximum 
cell  voltage  (MCV)  for  solid  oxide  fuel  cells  (SOFCs)  at  a  given  operating  condition.  These  fuels  make  a 
continuous  curve  in  the  ternary  coordinate  system.  A  fuel  map  is  established  by  developing  the  con¬ 
tinuous  fuel  curves  for  different  MCVs  at  the  same  operating  condition  and  representing  them  in  the 
carbon-hydrogen-oxygen  (C-H-O)  ternary  diagram.  Using  the  fuel  map,  the  effect  of  the  composition 
of  a  fuel  containing  carbon,  hydrogen,  oxygen,  and  inert  gas  atoms  on  the  MCV  of  SOFCs  can  be  easily 
studied.  In  addition  to  the  effect  of  the  fuel  composition,  the  graphical  representation  of  fuel  maps  can 
be  applied  to  study  the  effect  of  the  fuel  processors  on  the  MCV  of  SOFCs.  As  a  general  result,  among 
fuels  that  can  be  directly  utilized  in  SOFCs,  at  the  same  temperature  and  pressure,  the  one  located  at  the 
intersection  of  the  H-C  axis  and  the  carbon  deposition  boundary  (CDB)  curve  in  the  C-H-0  ternary  dia¬ 
gram,  provides  the  highest  MCV  for  SOFCs.  The  results  also  show  that  for  the  fuels  that  cannot  be  directly 
utilized  in  SOFC,  the  steam  reforming  fuel  processor  always  yields  a  higher  MCV  than  the  autothermal 
reforming  or  the  partial  oxidation  fuel  processors  at  the  same  inlet  fuel  temperature. 

©  2009  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

In  addition  to  the  high  electrical  efficiency,  the  fuel  flexibil¬ 
ity  is  one  of  the  most  important  advantages  of  solid  oxide  fuel 
cells  (SOFCs).  Various  types  of  fuels,  including  hydrogen,  carbon 
monoxide,  hydrocarbons,  alcohols  and  biomass  can  be  utilized  indi¬ 
rectly  or  even  directly  in  SOFCs.  Usually,  a  hydrogen-  and/or  carbon 
monoxide-rich  fuel  that  does  not  cause  any  carbon  deposition  over 
the  anode  catalyst  can  be  directly  utilized  in  SOFCs.  The  steam 
reforming  [1-3],  partial  oxidation  [4,5],  and  autothermal  reforming 
[6,7]  fuel  processors  are  widely  used  to  prepare  fuels  that  cannot 
be  directly  utilized  in  SOFCs.  In  addition  to  these  fuel  processors, 
anode  exit  gas  recirculation  [8,9]  can  be  applied  to  prevent  carbon 
deposition  over  the  anode  catalyst.  Although  a  broad  range  of  fuels 
and  fuel  processors  can  be  used  for  SOFCs,  their  effect  on  the  cell 
performance  can  vary  considerably.  In  this  paper,  the  effect  of  the 
composition  of  any  fuel  containing  carbon,  hydrogen,  oxygen,  and 
inert  gas  atoms  on  the  maximum  cell  voltage  (MCV)  of  SOFCs  is 
studied. 

Since  the  fuel  and  oxidant  are  utilized  along  the  anode  and  cath¬ 
ode  of  an  SOFC  (the  partial  pressures  of  fuel  and  oxidant  decrease 
along  the  anode  and  cathode  and  the  partial  pressure  of  prod¬ 
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ucts  increase  along  the  anode),  depending  on  the  cell  temperature 
distribution,  the  Nernst  voltage  and  internal  resistance  change  at 
different  places  of  a  cell.  In  most  cases,  especially  at  a  high  fuel 
utilization  ratio,  the  minimum  local  Nernst  voltage  and  the  max¬ 
imum  local  internal  resistance  take  place  at  the  end  of  the  anode 
(the  last  point  in  the  anode  along  the  fuel  flow  channel,  where  elec¬ 
trochemical  reactions  take  place).  In  this  study,  it  is  assumed  that 
the  fuel  utilization  ratio  is  sufficiently  high  so  that  the  minimum 
local  Nernst  voltage  and  the  maximum  local  internal  resistance  of  a 
cell  take  place  at  the  end  of  the  anode.  At  this  condition,  the  MCV  of 
ion-conducting  electrolyte  SOFCs  that  operate  at  a  specified  anode 
inlet  fuel  composition  and  temperature,  anode  outlet  fuel  temper¬ 
ature,  cell  operating  pressure,  oxygen  mole  fraction  at  the  cathode, 
and  fuel  utilization  ratio,  is  limited  to  the  minimum  local  Nernst 
voltage  and  can  be  thermodynamically  obtained  from  the  Nernst 
voltage  equation. 

Fuel  maps  are  developed  to  study  the  effect  of  the  inlet  fuel  com¬ 
position  on  the  MCV  of  SOFCs.  For  this  purpose,  a  thermodynamic 
model  is  developed  to  identify  the  fuels  that  would  yield  identi¬ 
cal  MCV  for  SOFCs  at  a  specified  operating  condition  of  anode  inlet 
and  outlet  fuel  temperatures,  cell  operating  pressure,  oxygen  mole 
fraction  at  the  cathode,  and  the  fuel  utilization  ratio.  These  fuels 
make  a  continuous  curve  in  the  ternary  coordinate  system.  A  fuel 
map  is  created  by  generating  the  continuous  fuel  curves  for  differ¬ 
ent  MCVs  at  the  same  operating  condition  and  representing  them  in 
the  carbon-hydrogen-oxygen  (C-H-O)  ternary  diagram.  It  is  nec¬ 
essary  to  draw  the  carbon  deposition  boundary  (CDB)  curve  in  a 
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Nomenclature 

C  carbon  atom  ratio 

C  consumption  rate  (mol  s-1 ) 

F  Faraday  constant  (=96485  Coulomb  mol-1 ) 

H  hydrogen  atom  ratio 

g  Gibbs  free  energy  (J  mol-1 ) 

h  enthalpy  (J  mol-1 ) 

i  electric  current  density  (A  M-2 ) 

I  electric  current  (A) 

Kp  equilibrium  constant 

MCV  maximum  cell  voltage  (V) 

h  molar  flow  rate  (mol  s-1 ) 

0  oxygen  atom  ratio 

p  pressure  (bar) 

P  production  rate  (mol  s-1 ) 

Ru  universal  gas  constant  (=8.314J  mol-1  K-1 ) 
T  temperature  (I<) 

Uf  fuel  utilization  ratio 

w  electric  power  (W) 

x  mole  fraction 

y  volumetric  mixing  ratio 

z  ratio  of  the  number  of  atoms 

Greek  letters 

a,p,y,8,£,\l/,co  defined  parameters 
A  difference 

Subscripts  and  superscripts 
0  standard  pressure  (=1  bar) 

B  Boudouard  reaction 

CD  carbon  decomposition  reaction 

CDB  carbon  deposition  boundary 

COR  carbon  monoxide  reduction  reaction 

i  inlet 

ig  inert  gas 

o  outlet 

SR  steam  reforming  reaction 

WGS  water  gas  shift  reaction 


fuel  map  to  determine  the  region  of  a  fuel.  The  region  below  the 
CDB  curve  belongs  to  fuels  that  can  be  directly  utilized  in  SOFCs.  In 
this  paper,  this  region  is  called  “direct  fuel  utilization  region”.  The 
region  above  the  CDB  curve  belongs  to  fuels  that  should  be  pro¬ 
cessed  before  being  utilized  in  SOFCs  to  prevent  carbon  deposition 
over  the  anode  catalyst.  In  this  paper,  this  region  is  called  “indirect 
fuel  utilization  region”.  After  finding  the  CDB  curve,  the  fuel  map 
can  be  easily  applied  to  study  the  effect  of  the  composition  of  any 
fuel  containing  carbon,  hydrogen,  oxygen  and  inert  gas  atoms  on 
the  MCV  of  SOFCs. 

2.  Model  development 

At  the  operating  temperatures  and  pressures  of  SOFCs,  the 
species  of  primary  importance  of  the  equilibrium  products  of  a 
fuel,  located  in  the  C-FI-0  ternary  diagram,  are  carbon,  C,  in  solid 
phase  and  a  mixture  of  H2,  CO,  CFI4,  FI20  and  C02  in  the  gaseous 
phase.  The  species  of  secondary  importance,  which  are  present  in 
minor  amounts,  are  C2FI6,  C2FI4,  and  C2FI2.  Species  such  as  CFI3OFI, 
HCHO,  C2H5  OH,  and  Ci0H22  are  present  in  amounts  of  several  orders 
of  magnitude  smaller  than  the  species  of  secondary  importance 
[10,11].  The  chemical  species  of  the  equilibrium  products  of  300 
different  fuels,  located  in  the  C-H-0  ternary  diagram,  including 
alkanes,  alcohols,  alkenes,  alicyclic  hydrocarbons,  dimethyl  ether, 


biogas  and  coke  oven  gas  have  also  been  calculated  by  Sasaki  et  al., 
via  the  Gibbs  free  energy  minimization  method  [12].  Their  results 
for  the  major  and  minor  constituents  of  equilibrium  products  of 
fuels  have  a  good  agreement  with  the  results  presented  in  refer¬ 
ences  [10,11]. 

In  the  present  model,  all  fuels  located  in  the  C-H-0  ternary  dia¬ 
gram  are  assumed  to  be  in  thermodynamic  equilibrium  and  only 
the  important  species  of  H2,  CO,  CH4,  H20  and  C02  are  taken  into 
account  for  the  gaseous  phase  of  equilibrium  products  of  the  fuels. 
However,  inert  gases  such  as  nitrogen  may  be  present  in  this  phase. 
Neglecting  the  amounts  of  the  minor  constituents  in  the  model 
does  not  affect  the  results  within  the  accuracy  of  the  equilibrium 
constant  data  [11]. 

Before  presenting  the  model,  the  mole  fraction  of  species  in  the 
gaseous  phase  of  equilibrium  product  of  a  fuel  (hereafter  called  just 
fuel)  is  obtained  by  dividing  the  partial  pressures  of  H2,  CO,  CH4, 
H20,  and  C02  by  the  difference  of  the  total  gas  pressure  and  the 
inert  gases  partial  pressure.  Using  the  Dalton’s  law: 


*h2  +* *co  +  *ch4  +  *h2o  +  *C02  =  1  (1) 

The  location  of  a  gaseous  fuel  in  the  C-H-0  ternary  diagram  can  be 
determined  from  [11]: 


*co2  +  *co  +  *ch4 

2xh2  +  3xCo2  +  2xCo  +  3xH2o  +  5xCh4 

2xh2  +  2xh2o  +  4xCh4 
2xh2  +  3xCo2  +  2xCo  +  3xh2o  +  5xCh4 

2xCo2  +  *co  +  *h2o 

2xh2  +  3xCq2  +  2xCo  +  3xh2o  +  5xCh4 


(2) 

(3) 

(4) 


2.1.  Carbon  deposition  boundary 


Depending  on  the  fuel  composition,  temperature  and  pressure, 
the  presence  of  the  solid  carbon  in  a  gaseous  mixture  is  possible. 
The  solid  carbon  may  form  on  the  anode  catalyst  and  deactivate 
it  for  electrochemical  reactions  [13,14].  As  a  preliminary  step  to 
develop  fuel  maps,  the  CDB  curve  should  be  determined  based 
on  thermodynamic  equilibrium  or  considering  the  influence  of 
the  anode  catalyst  using  a  detailed  kinetics  modeling  or  experi¬ 
ment.  These  methods  have  been  studied  extensively  by  researchers 
[11,12,15-24].  To  obtain  the  CDB  curve  for  fuel  maps  in  this  paper, 
the  thermodynamic  equilibrium  method  suggested  by  Broers  et  al. 
[11]  was  simplified. 

In  thermodynamic  equilibrium,  the  CDB  curve  can  be  deter¬ 
mined  considering  that  the  solid  carbon  can  be  formed  by  the  three 
reactions  of  carbon  decomposition  (Rl),  CO  reduction  (R2)  and  the 
Boudouard  reaction  (R3)  [25,26].  On  the  CDB  curve,  the  net  amount 
of  solid  carbon  production  and  consumption  by  these  reactions  is 
zero  [27]. 

CH4'^C  +  2H2  (Rl ) 


Kp.com 


XH2,CDb(P  Pig) 
*ch4,cdb 


h2  +  coV^c  +  h2o 


(5) 

(R2) 


Rp,cor  (T)  = 


*h2o,cdb 

*H2,CDB*CO,CDb(P  -  Pig) 


(6) 


Kp  B 

2C0^>C  +  C02 


(R3) 


I<p,b(T) 


*co2,cdb 

XCO,CDB (P  ~  Pig) 


(7) 


The  mole  fractions  of  CH4,  H20  and  C02  at  the  given  fuel  tem¬ 
perature,  total  pressure,  and  partial  pressure  of  the  inert  gases  are 
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determined  from  Eqs.  (5)-(7),  as  follows: 


XH2,CDb(P  Pig) 

Xch4,cdb  =  KpCD{T) 

(8) 

^H20,CDB  =  ^p,COr(^)^H2,CDB^CO,CDb(P  -  Pig) 

(9) 

*C02,CDB  —  ^p,b(^)^C0,CDb(P  —  Pig) 

(10) 

By  substituting  Eqs.  (8)-(10)  into  Eq.  (1),  the  mole  fraction  of  H2 
can  be  determined  as  a  function  of  the  mole  fraction  of  CO. 


*h2,cdb  =  0.5I<pcd{T) 


+  Kp,coR(r)Xco,( 


■0 


-*CO,CDB  js 

(P-Pig)  “  KP’b 


vr  r  (t)xco,cdb) 

1  + - ^ , 

^  Kp,CD  (^)((p-pig)  +^p,cor(7’)^co,cdb^ 


-1 


(11) 


The  composition  of  fuels  that  makes  the  CDB  curve  in  the  C-H-0 
ternary  diagram  is  determined  using  Eqs.  (8)— (11 )  by  changing  the 
mole  fraction  of  carbon  monoxide  from  0  to  Xco,CDB,max- 


o  <  *CO,CDB  <  ^CO, CDB, max  = 

2.2.  Maximum  cell  voltage 


y'l  +4Kp,B(T)(p-pig)-  1 
2KP,B(T)(P-Pig) 


(12) 


To  obtain  the  MCV  of  SOFCs,  it  is  assumed  that  the  fuel  utilization 
ratio  is  sufficiently  high  so  that  the  minimum  local  Nernst  voltage 
and  the  maximum  local  internal  resistance  in  a  cell  take  place  at 
the  end  of  the  anode.  From  the  Nernst  voltage  equation: 

—RuTq  , 


MCV  = 


2  F 


-[oi  +  ln(e)) 


fS,o(7'o)-fS7(7’o)-0.5i°  (T0) 


RuTo 


-0.5  In (x0  p) 


£  = 


*H20,o 

*H2,o 


(13) 

(14) 

(15) 


To  determine  the  MCV,  Eq.  (15)  should  be  expanded.  For  this 
purpose,  the  molar  flow  rate  of  each  species  in  the  outlet  fuel  is 
derived  as  follows: 


^CH4,o  =  ^CH4,i  -  Qh4 

nc02,o  =  Hco2,i  +^co2 
nco,o  =  nC0)i  +  Cch4  -  Pc o2 
^H2,o  =  tlu 2,i  +  3Cch4  +  Pco2  -  ch2 
^H20,o  =  ”H20,i  -  Qh4  -  Pc02  +  Ch2 


(16) 

(17) 

(18) 

(19) 

(20) 
(21) 


CH4  +  H20<£5c0  +  3H2 


Kp,sr  (T)  = 


*co*h2(P  -  Pig) 


CO  +  H20  #^sco2  +  H2 

Kp,wcs(T)=^f± 

*co*h2o 


(R5) 

(22) 

(R6) 

(23) 


The  difference  of  the  molar  flow  rate  of  the  outlet  fuel  and  the  inert 
gases  is  obtained  from  Eqs.  (16)-(20)  as  follows: 


(ii0  -  hig  „)  =  (hi  -  hig  0  +  2CCh4 


(24) 


By  dividing  Eqs.  (16)-(20)  by  Eq.  (24),  the  mole  fractions  of  CH4, 
C02,  CO,  H2,  and  H20  in  the  outlet  fuel  are  obtained. 


*C02,o  = 


*CO,o  = 


*CH4,i  CcH4/(hj  ^ig,i) 

1  +  2CcH4/(hi  —  h\ g  j) 

*C02,i  +l)C02/(^i  —  ^ig,i) 

1  +2CCH4/(hi  -fiigJ) 

*co,i  +  CCH4/(ni  -  h[g  i)  -  Pco2/(hi  -  nig>i) 
1  +  2CCh 4/{jii  ~  tiigj) 


(25) 

(26) 

(27) 


*H2,o  — 


*H2,i  +  3CCH4/(hi  -  higii)  +  Pco2/(hi  -  higji)-  CH2/(hi  -  hig>i) 


1  “f  2CcH4/(^i  ^ig,i 


(28) 


*H20,o  = 


*H20,i  CcH4/(hj  ^ig,i)  Pc02/{tl[  ^igj) +Q-I2 /(^i 


1  T  2 Cq-i4  / ( hi  hig  i ) 


(29) 


By  substituting  Eqs.  (28)  and  (29)  into  Eq.  (15),  this  equation  is 
expanded  as  follows: 


£  = 


*H2o,i  ~  CCH4/(iii  -  higii)  -  PCo2/(hj  -  higii)  +  CH2/(hj  -  hig, j) 
*H2,i  T  3CcH4/(fIi  —  ^ig,i)  "f  Pco2/(^i  —  ^ig,i)  —  1-H2  /(^i  —  ^ig,i) 


(30) 

The  next  step  is  to  derive,  Ccu4/{h[  -  hig  i),  C\\2/{hx  -  hjgj)  and 
Pco2/(«i  -nig,i). 

By  utilizing  the  hydrogen  along  the  anode  and  increasing  the 
temperature  of  fuel,  reaction  (R5)  is  shifted  to  the  right-side  to  con¬ 
sume  almost  all  methane  content  of  the  outlet  fuel  from  the  anode; 
thus,  from  Eq.  (25): 


Qh4 


=  *CH. 


M  ,lig,i 


4d 


(31) 


Cu2/{h[  -  hig ?i)  can  be  obtained  from  the  fuel  utilization  ratio  equa¬ 
tion,  Eq.  (32)’ [28-30],  as  Eq.  (33). 


Uf  = 


[(^H2,i  -  ^h2,o)(^h2o( Ji)  -  hH2(W_-  0.5ho2(Ti))  +  (hCo,i  7  ”co,o)(^co2 (Tj)  -  ^co(^i)  -  0.5hO2(ri)) 

+  (^CH4,i  -  rcCH4,o)(frc02(^i)  +  2hH2o(^i)  ~  ^CH4(^i)  -  2ho2(Ti))] 

[hu2,iChH2o(Ti)  -  hH2(Ti)  -  0.5ho2_(Ti))  +  nco,i(hco2 (Ti )  -  hcoffi)  -  0.5hO2(Ti)) 

+^CH4,i(^C02(l'i)  +  2ftH2o(’Ti)  -  ^CH4(7i)  -  2ho2(^i))] 


(32) 


where  Ch2  and  CCh4  are,  respectively,  the  rates  of  the  hydrogen  and 
methane  consumption  by  reactions  (R4)  and  (R5)  and  Pc o2  is  the 
rate  of  the  carbon  dioxide  production  by  reaction  (R6). 

H2  +  02-^  H20  +  2e" 


(hi  -Hhigii) = Uf  + Vco:' + f  ^ 


(R4) 


8-y 


(33) 
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where, 

hH2o(7'i)-hH2(Ti)-0.5ho2(ri) 
P  RuTo 

hco2(7'i)-hco(ri)-0.5ho2(Ti) 


Y  = 


<5  = 


RuTo 

hcofJiH  2hH2o(Ti)  ~  haulTj)-  2/^(7)) 

RuTo 


(34) 

(35) 

(36) 


By  solving  Eqs.  (1),  (22)  and  (23)  for  the  inlet  fuel  to  the  anode, 
the  mole  fraction  of  CO2,  CO,  and  CH4  are  determined  as  follows: 

Kp,wcs(ri)x2  .(1  -  XH2o,i  -  XH2,i) 

XC02,i  =  - : - - 72 

2  XH  /P-PigJ 

^p,WGs(^i)^H2o,i  +  xH2,ixH20,i  H  Kp>SR(7j) 


(41) 


xCO,i  = 


xH2,ixH20,i(l xH20,i xH2,i) 


/Cr 


rp,WGs(l’i)^H20,i  +  XH2xH20,i  +  (xh2j(P  “  Pig)  /Kp,SR.{Ti)) 


(42) 


To  obtain  PCo2 /(hi  -  hig  i),Eq.  (23)  is  developed  for  the  outlet  fuel 
from  the  anode.  Depending  on  the  equilibrium  constant  of  the  water 
gas  shift  reaction  at  the  temperature  of  the  outlet  fuel,  Pc o2/(hi  - 
hig  j)  is  determined  from  the  following  equation: 


XH2?i(P Pig)2/^P,SR(^i)(^  -xH20,i  -XH2,i) 


CH4,‘  Rp,WCs(Ti)X^0  +XH2,iXH2o,i  +  (*h2  j(P -Pig)2/^p,SR(f'i)) 

(43) 


PcOy 


("i  -  tlig,i) 


(*CO,i  +  XcH4.i)  (  Xch4 , i  (  —jf~  ~  2  j  -Xn20,i  ~  Of  fxH2J  +  ^Xco  i  +  -gXCH4,i 


+  *C02,i  (  *H2,i  +xCH4,i  (  j  ~  Of  (  XH2,i  +  ^XCOjj  +  ^XCH4j 


x  (xC02,i  +  xCO,i  +  xCH4,i)  +  1 
'  2 


/  2xCH4,i+  (jj  -  l) 


'co2 


(hi  -  hig  i)  2 


+  'It, 


Re o2  ifr  \ 

(hi  -  hig,i)  2  +  V  (  2  1  ® 


(fCp,wcs(r0)  =  i) 
(iCp,wcs(r0)  >  i) 
(KP,wcs(r0)  <  i) 


(37) 


where, 


l/f  =  P 

y 


xH2,i  +  xC02,i  +  xCH4,i  ~  Uf  (  xH2,i  +  j^xCO,i  +  ^xCH4,ij  +  (  ^  ~  3  )  xCH4,i 


2xCH4,i  +  -  Xc°2  4  +^P,WGs(lo)  ^1+^-1^  (xCO,i  +  xCH4,i)^ 


1  -  ^p,WGs(^o) 


P 


y(i-^P,wGs(r0)) 


Kp , WGS ( lo )(xCO , i  +  xCH4,i)  ^xCH4,i  ^  ^  ~  2^  -  xH20,i  “  ^xH2,i  +  ^XCO,i  +  ^xCH4,i^  ^ 


+xC02,i  (  xH2,i  +  xCH4,i  (  ^  J  ~  Uf  fxH2,i  +  ^XCO,i  +  ^xCH4,i 


(38) 


(39) 


Therefore,  at  a  given  operating  condition  of  anode  inlet  and  out¬ 
let  fuel  temperatures;  cell  operating  pressure;  oxygen  mole  fraction 
at  the  cathode;  and  fuel  utilization  ratio,  the  MCV  depends  only  on 
the  anode  inlet  fuel  composition. 


2.3.  Developing  fuel  map 


In  this  step,  the  equations  required  to  find  the  inlet  fuels  to  the 
anode  that  yield  identical  MCV  at  the  given  operating  condition  are 
determined. 

At  a  given  operating  condition  and  a  specific  MCV,  the  mole 
fraction  of  H20  in  the  inlet  fuel  can  be  found  from  Eq.  (13): 


xH20,i  =  (  xH2,i  +3xCH4,i  +  ^ 

(  —RuTo 


Pc  02 


Ch2 


77;  -  77 


ig,i 


x  exp 


V  2F  MCV 


-  xH2,i  -  2xCH4,i 


(40) 


By  changing  the  mole  fraction  of  H2  in  the  inlet  fuel  from  0  to  1, 
the  compositions  of  the  inlet  fuels  that  yield  an  identical  MCV  for 
SOFCs  at  the  given  operating  condition  are  determined  from  Eqs. 
(40)-(43).  Among  these  fuels,  those  located  below  the  CDB  curve 
in  the  C-H-0  ternary  diagram  are  acceptable.  These  fuels  create  a 
continuous  curve  in  the  ternary  coordinate  system.  A  fuel  map  is 
established  by  developing  the  continuous  fuel  curves  for  different 
MCVs  at  the  same  operating  condition  and  representing  them  in 
the  C-H-0  ternary  diagram. 

3.  Results  and  discussion 

The  locations  of  some  typical  fuels,  presented  in  Tables  1  and  2, 
are  shown  in  Fig.  1.  Most  of  these  fuels  are  above  the  carbon  depo¬ 
sition  curve  in  the  C-H-0  ternary  diagram  (in  the  indirect  fuel 
utilization  region)  and  cannot  be  directly  utilized  in  SOFC.  Indeed, 
they  should  be  transferred  to  below  the  CDB  curve  by  combining 
them  with  a  mixture  (or  a  substance)  located  below  the  CDB  curve 


Table  1 

A  typical  composition  of  some  gaseous  fuels. 


S.  Farhad,  F.  Hamdullahpur  /  Journal  of  Power  Sources  191  (2009)  407-416 


411 


Compound  mol% 


NGa  [31] 

ADGb  [32] 

SRMC  [33] 

POMd  [34] 

FBCGe  [35] 

BPGf  [36] 

ch4 

88.5 

60.8 

2.4 

1.0 

4.6 

4 

c2h6 

4.8 

- 

- 

- 

- 

- 

c3h8 

1.6 

- 

- 

- 

- 

- 

c4h10 

0.7 

- 

- 

- 

- 

- 

h2 

- 

- 

52.9 

35.8 

28.0 

17 

CO 

- 

- 

0.50 

18.9 

33.0 

13 

h2o 

- 

0.01 

30.4 

3.7 

16.8 

15 

co2 

- 

34.8 

13.0 

0.8 

15.0 

11 

n2 

4.4 

2.4 

0.8 

39.8 

0.6 

40 

02 

- 

1.5 

- 

- 

- 

- 

a  Natural  gas. 
b  Anaerobic  digestion  gas. 

c  Gas  produced  by  the  steam  reforming  of  methane. 
d  Gas  produced  by  the  partial  oxidation  of  methane. 
e  Gas  produced  by  the  fluidized  bed  coal  gasifier. 
f  Biomass  produced  gas. 


(e.g.  with  steam  for  the  case  of  the  steam  reforming  or  with  oxygen 
for  the  case  of  the  partial  oxidation  fuel  processors).  In  the  C-H-0 
ternary  diagram,  any  possible  combination  of  two  different  mix¬ 
tures  (or  substances)  lies  on  a  line  that  connects  the  location  of 
these  two  mixtures.  Depending  on  the  location  of  the  new  mixture 
on  this  line,  the  volumetric  mixing  ratio  of  mixture  2  to  mixture  1 
is  obtained  from  Eq.  (44). 

where  Z\\2  is  the  ratio  of  the  number  of  atoms  in  mixture  1  to  mix¬ 
ture  2,  and  Qi2  is  the  carbon  atom  ratio  of  the  new  mixture.  Hence, 
the  volumetric  mixing  ratio  of  a  fuel  located  above  the  CDB  curve 
and  a  mixture  (or  a  substance)  below  this  curve,  which  makes  a  fuel 
that  can  be  utilized  in  the  anode,  can  be  determined  from  Eq.  (44). 

3.1  Carbon  deposition  boundary 

In  Fig.  2,  the  effect  of  the  fuel  temperature  and  pressure  on  the 
CDB  curve  is  shown.  According  to  Fig.  2a,  increasing  the  temper- 


(a) 


O 


ATOM  PERCENT  O 

Fig.  1.  The  location  of  some  fuels  in  the  C-H-0  ternary  diagram. 


Fig.  2.  The  CDB  curve  in  the  C-H-0  ternary  diagram:  (a)  effect  of  temperature  and 
(b)  effect  of  pressure. 
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Fig.  3.  Fuel  map  for  the  MCV  of  SOFCs  operating  at  p  =  1  atm,  Uf  =  0.85,  A T=  200  K,  x02  =  0.21 ,  pig  =  0  and  (a)  T{  =  700  K,  (b)  Tj  =  800  K,  (c)  f  =  900  K,  (d)  f  =  1000  K,  (e)  f  =  1100  K, 
and  (f)Ti  =  1200 IC 
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Table  2 

A  typical  composition  of  some  solid  fuels  [37]. 


Feedstock 

wt%  (dry) 

C 

H 

0 

N 

S 

Ash 

MSWa 

37.7 

5.1 

30.1 

1.2 

0.1 

25.8 

Waste  wood 

41.6 

4.9 

36.2 

0.3 

0.01 

17 

Scrap  tire 

71.5 

6.9 

0.7 

0.6 

1 

19.3 

Peat 

50.7 

5.3 

34.7 

1.3 

1.3 

6.7 

Lignite 

60 

6 

20 

2 

2 

10 

Bituminous  coal 

73.9 

4.9 

9.2 

1.4 

0.8 

9.8 

Anthracite 

89.5 

3.1 

2.1 

1.1 

1.1 

3.1 

a  MSW:  municipal  solid  waste. 


Table  3 

The  fuels  that  yield  the  MCV  =  0.82  V  at  7)  =  900  K,  A T=  200  K,  p  =  1  atm,  Uf  =  0.85,  and 
Xq2  =  0.21. 


ature  causes  the  direct  fuel  utilization  region  to  enlarge  on  the 
right-  and  shrink  on  the  left-side  of  the  C-H-0  ternary  diagram. 
At  temperatures  of  more  than  1200  K,  the  CDB  curve  is  almost  a 
straight  line  joining  the  points  H2  and  CO  together. 

As  shown  in  Fig.  2b,  the  effect  of  decreasing  the  pressure  on  the 
CDB  curve  is  similar  to  the  effect  of  reducing  temperature.  In  fact,  by 
increasing  the  pressure,  the  direct  fuel  utilization  region  shrinks  on 
the  right-  and  enlarges  on  the  left-side.  Based  on  Eqs.  (8)— (11 ),  the 
effect  of  inert  gases  on  the  CDB  curve  is  equivalent  to  decreasing 
the  total  pressure  of  the  gaseous  fuel.  This  means  that  a  gaseous 
fuel  containing  20  mol%  of  inert  gases  is  equivalent  to  the  same  fuel 
without  any  inert  gas  and  a  total  pressure  of  20%  less  than  the  fuel 
with  inert  gas. 
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Fig.  5.  Effect  of  the  inlet  fuel  temperature  on  the  highest  MCV  at  different  fuel 
utilization  ratios  (p  =  1  atm,  AT  =  200  K,  xq2  =  0.21 ). 


3.2.  Maximum  cell  voltage 

The  predictions  obtained  from  the  current  model  for  some  fuels 
that  yield  the  same  MCV  are  shown  in  Table  3  for  the  inlet  and 
outlet  fuel  temperatures  of  Tx  =  900  K  and  T0  =  1100  K,  respectively 
(AT  =  r0  -  T[  =  200KJ;  the  cell  operating  pressure  ofp  =  latm;  the 
fuel  utilization  ratio  of  Df  =  0.85;  and  the  oxygen  mole  fraction  of 
Xo2  =  0.21.  According  to  this  table,  the  presence  of  methane  in  the 
inlet  fuel  does  not  necessarily  lead  to  an  increase  in  the  MCV. 

In  Fig.  3,  the  fuel  maps  to  predict  the  effect  of  fuel  composi¬ 
tion  on  the  MCV  are  shown  for  different  inlet  fuel  temperatures  of 
700  K,  800  K,  900  K,  1000  K,  1100  K,  and  1200  K  atp  =  1  atm,  Uf  =  0.85, 
AT =200  K,  Xq2  =  0.21 ,  and  pig  =  0.  It  can  be  seen  from  this  figure  the 


3 


4 

2 


T(K) 


T(K) 


Fig.  4.  Effect  of  the  temperature  on  the  mole  fraction  of  methane  in  the  fuel  that 
yields  the  highest  MCV  for  SOFCs  at  different  cell  operating  pressures. 


Fig.  6.  Effect  of  the  inlet  fuel  temperature  on  the  highest  MCV  at  different  fuel 
temperature  rises  along  the  anode  (p  =  1  atm,  Uf  =  0.85,  Xq2  —  0.21 ). 
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Fig.  7.  Effect  of  the  (a)  temperature  of  the  inlet  fuel  to  the  anode  at  MCV  =  0.76  V,  p  =  1  atm,  Uf  =  0.85,  A T=  200  K,  Xq2  =  0.21 ,  pig  =  0,  (b)  cell  operating  pressure  at  MCV  =  0.82  V, 
Tj  =  900  K,  Uf  =  0.85,  AT  =  200  K,  Xo2  =  0.21 ,  pig  =  0,  (c)  fuel  temperature  rise  along  the  anode  at  MCV  =  0.82  V,  Tj  =  900  K,  p  =  1  atm,  Uf  =  0.85,  Xo2  =  0.21 ,  pig  =  0,  (d)  fuel  utilization 
ratio  at  MCV  =  0.82  V,  f  =  900  K,  p  =  1  atm,  AT  =  200  K,  xo2  =  0.21 ,  pig  =  0,  (e)  oxygen  mole  fraction  in  cathode  at  MCV  =  0.82  V,  T,  =  900  K,  p  =  1  atm,  Uf  =  0.85,  A T=  200  K,  pig  =  0, 
and  (f)  nitrogen  mole  fraction  (as  an  inert  gas)  at  MCV  =  0.86  V,  7)  =  800  K,  p  =  1  atm,  Uf  =  0.85,  A T=  200  K,  xq2  =  0.21 ,  on  the  continuous  fuel  curves. 
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type  of  fuel  has  a  significant  effect  on  the  MCV.  The  following  gen¬ 
eral  results  can  be  obtained  for  the  fuels  in  the  direct  and  indirect 
fuel  utilization  regions  in  the  C-H-0  ternary  diagram. 

3.2  A.  Direct  fuel  utilization  region 

Among  the  fuels  with  the  same  percentage  of  carbon  atom,  the 
one  with  a  higher  percentage  of  hydrogen  atom  yields  a  higher  MCV 
at  a  given  operating  condition.  If  the  location  of  a  fuel  is  either  below 
or  at  right  of  another  fuel  in  the  C-H-0  ternary  diagram,  it  yields  a 
lower  MCV. 

Among  the  fuels  that  are  at  the  same  temperature  and  pressure, 
the  highest  MCV  is  yielded  by  the  one  located  at  the  intersection  of 
the  CDB  curve  and  the  H-C  axis  in  the  C-H-0  ternary  diagram.  This 
fuel  contains  only  H2  and  CH4;  however,  a  small  quantity  of  water 
should  be  added  to  start  the  electrochemical  reaction  in  anode.  For 
this  fuel,  the  mole  fraction  of  CH4  as  a  function  of  the  inlet  fuel  tem¬ 
perature  at  different  cell  operating  pressures  is  shown  in  Fig.  4.  It  is 
observed  that  decreasing  the  inlet  fuel  temperature  and  increasing 
the  cell  operating  pressure  leads  to  an  increase  in  the  mole  fraction 
of  CH4.  At  inlet  fuel  temperatures  of  more  than  1200 1<  and  a  cell 
operating  pressure  of  the  order  of  1  atm,  almost  pure  hydrogen  is 
the  fuel  that  yields  the  highest  MCV  for  SOFCs.  The  effect  of  the  inlet 
fuel  temperature  on  the  highest  MCV  at  different  fuel  utilization 
ratios  and  fuel  temperature  rises  along  the  anode  (T0  -  7/)  are  illus¬ 
trated  in  Figs.  5  and  6,  respectively.  According  to  these  figures,  for 
each  100  K  increase  in  the  inlet  fuel  temperature,  the  highest  MCV 
reduces  around  0.04  V.  Increasing  the  fuel  utilization  ratio  and  fuel 
temperature  rise  along  the  anode  causes  a  decrease  in  the  highest 
MCV. 

In  many  cases,  for  a  fuel  in  the  direct  fuel  utilization  region, 
the  MCV  can  be  increased  by  mixing  it  with  a  fuel  which  is  in  the 
indirect  fuel  utilization  region.  For  example,  if  a  mixture  of  85  mol% 
biomass  produced  gas  (BPG)  and  15  mol%  natural  gas  (NG),  see  Fig.  1 , 
is  utilized  in  an  SOFC  instead  of  utilizing  only  the  BPG,  based  on  the 
fuel  map  presented  in  Fig.  3e,  the  MCV  increases  from  0.73  V  to 
0.75  V  at  the  operating  condition  of  Tx  =  1100  K,  p  =  1  atm,  Uf  =  0.85, 
AT  =  200  K,  andxo2  =0.21. 

3.2.2.  Indirect  fuel  utilization  region 

In  the  indirect  fuel  utilization  region,  the  fuel  composition  and 
the  type  of  the  fuel  processor  affects  the  MCV  of  SOFCs.  Among 
the  fuels  with  the  same  percentage  of  carbon  atom,  a  fuel  with  a 
higher  percentage  of  hydrogen  atom  yields  a  higher  MCV  at  a  given 
operating  condition,  regardless  of  the  type  of  the  fuel  processor. 

We  decided  to  use  an  example  to  explain  the  effect  of  different 
types  of  fuel  processors  on  MCV.  In  this  example,  CH4  at  a  tem¬ 
perature  of  900  K  is  considered  as  the  inlet  fuel  to  the  anode.  Using 
Fig.  2a,  we  can  determine  that  this  fuel  is  in  the  indirect  fuel  utiliza¬ 
tion  region  and  should  be  processed  before  it  can  be  utilized  in  the 
anode.  If  partial  oxidation  is  selected  as  a  fuel  processor,  according 
to  Eq.  (44)  and  Fig.  3c,  at  least  0.86  moles  of  02  should  be  added  to 
each  mole  of  CH4  to  cross  the  CDB  curve  which  will  then  result  in  an 
MCV  value  around  0.83  V.  If  steam  reforming  is  selected  as  the  fuel 
processor,  at  least  1.4 moles  of  H20  should  be  added  to  each  mole 
of  CH4  and  the  MCV  will  be  around  0.84  V.  For  the  auto-thermal 
reforming  fuel  processor,  depending  on  how  much  H20  and  02  are 
added  to  CH4,  the  MCV  varies  from  0.83  V  to  0.84  V.  The  dry  reform¬ 
ing  of  CH4  at  900 1<  is  possible  by  adding  almost  6.7  moles  of  C02  to 
each  mole  of  CH4;  however,  the  MCV  reaches  around  0.8  V.  If  four 
moles  of  H2  are  added  to  each  mole  of  CH4  at  900  K,  the  MCV  cor¬ 
responding  to  the  new  mixture  will  be  around  0.85  V.  The  effect  of 
the  anode  exit  gas  recirculation  on  the  MCV  is  the  same  as  the  par¬ 
tial  oxidation  with  oxygen.  As  a  general  result,  for  any  fuel  in  the 
indirect  fuel  utilization  region,  the  steam  reforming  fuel  processor 
yields  a  higher  MCV  than  the  partial  oxidation  and  auto-thermal 
reforming  fuel  processors  at  a  specific  operating  condition. 


3.3.  Effect  of  the  operating  condition 

To  predict  the  effect  of  the  operating  condition  of  a  cell  on  the 
location  of  the  continuous  fuel  curves  in  fuel  maps,  a  parametric 
study  is  performed  and  the  results  are  shown  in  Fig.  7.  According 
to  Fig.  7a,  decreasing  the  temperature  of  the  inlet  fuel  to  the  anode 
shifts  the  continuous  fuel  curves  to  the  right  of  the  fuel  map.  This 
means  that  fuels  with  a  lower  hydrogen  atom  percentage  (with  the 
same  percentage  of  carbon  atom)  can  still  yield  a  high  MCV  if  the 
inlet  fuel  temperature  is  reduced.  Based  on  Fig.  7b,  by  decreasing  the 
cell  operating  pressure,  the  continuous  fuel  curves  tend  to  move  to 
the  left.  As  observed  in  Fig.  7c  and  d,  by  decreasing  the  fuel  temper¬ 
ature  rise  along  the  anode  and  the  fuel  utilization  ratio,  the  contin¬ 
uous  fuel  curves  shift  to  the  right  of  the  fuel  map.  By  decreasing  the 
oxygen  mole  fraction  in  cathode,  the  continuous  fuel  curves  shift  to 
the  left  of  the  fuel  map  (Fig.  7e).  According  to  Fig.  7f,  there  is  no  sig¬ 
nificant  change  in  the  location  of  continuous  fuel  curves  between 
a  fuel  without  any  inert  gas  and  a  fuel  with  75  mol%  inert  gas. 

4.  Conclusions 

At  any  operating  condition,  a  fuel  map  can  be  developed  for 
SOFCs  to  predict  the  effect  of  the  composition  of  any  fuel  containing 
carbon,  hydrogen,  oxygen,  and  inert  gas  atoms  on  the  MCV.  The  fuel 
maps  are  also  useful  for  selecting  an  appropriate  fuel  processor  for 
the  fuels  that  cannot  be  directly  utilized  in  SOFCs  due  to  the  carbon 
deposition  problem.  Some  of  the  general  results  obtained  in  this 
study  are  as  follows: 

•  The  composition  of  a  fuel  has  a  significant  effect  on  the  MCV  of 
SOFCs. 

•  Among  fuels  with  the  same  percentage  of  carbon  atom,  a  fuel  with 
a  higher  percentage  of  hydrogen  atom  yields  a  higher  MCV  at  a 
given  operating  condition. 

•  In  the  direct  fuel  utilization  region,  if  the  location  of  a  fuel  is  either 
below  or  right-side  of  another  fuel  in  the  fuel  map,  it  yields  a  lower 
MCV. 

•  Among  the  fuels  in  the  direct  fuel  utilization  region,  at  the  same 
temperature  and  pressure,  the  highest  MCV  is  always  obtained 
from  a  fuel  located  at  the  intersection  of  the  CDB  curve  and  the 
H-C  axis  in  the  fuel  map. 

•  In  the  indirect  fuel  utilization  region,  the  composition  of  a  fuel 
and  the  type  of  the  fuel  processor  affect  the  MCV. 

•  For  any  fuel  in  the  indirect  fuel  utilization  region,  the  steam 
reforming  fuel  processor  yields  a  higher  MCV  than  the  partial 
oxidation  and  auto-thermal  reforming  fuel  processors  at  a  given 
operating  condition  of  a  cell. 

•  Decreasing  the  cell  operating  pressure  and  the  oxygen  mole  frac¬ 
tion  (in  cathode)  shifts  continuous  fuel  curves  to  the  left  of  the 
fuel  map. 

•  Decreasing  the  temperature  of  the  inlet  fuel  to  the  anode;  the  fuel 
temperature  rise  along  the  anode;  or  the  fuel  utilization  ratio, 
shifts  the  continuous  fuel  curves  to  the  right  of  the  fuel  map. 

•  The  partial  pressure  of  inert  gases  does  not  have  any  significant 
effect  on  the  location  of  continuous  fuel  curves. 
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